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Fundamental parameters of RR Lyrae stars from 
multicolour photometry and Kurucz atmospheric models. 
I. Theory and practical implementation 
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ABSTRACT 

A photometric calibration of Kurucz static model atmospheres is used to obtain the 
following parameters of RR Lyrae stars: variation of stellar angular radius ^, effective 
temperature To, gravity g^ as a function of phase, interstellar reddening E{B — V) 
towards the star and atmospheric metallicity M . Photometric and hydrodynamic con- 
ditions are given to find the phases of pulsation when the quasi-static atmosphere 
approximation (QSAA) can be applied. The QSAA is generalized to a non-uniformly 
moving spherical atmosphere, and the distance d, mass A4 and atmospheric motion 
are derived from the laws of mass and momentum conservation. To demonstrate the 
efficiency of the method, the UBV(RI)c photometry of SU Dra was used to de- 
rive the following parameters: [M] = -1.60 ± .10 dcx, E{B - V) ^ 0.015 ± .010, 
d = 663 ± 67 pc, M = (0.68 ± .O3)A^0, equilibrium luminosity Loq = 45.9 ± 9.3Lq 
and Tcq == 6813 ± 20 K. 

Key words: hydrodynamics - stars: atmospheres - stars: fundamental parameters - 
stars: individual: SU Dra - stars: oscillations - stars: variables: RR Lyrae. 



1 INTRODUCTION 

Although the determination of the fundamental parameters 
of RR Lyrae (RR) stars is interesting in itself, it is also im- 
portant from a practical point of view because RR stars play 
a considerable role in establishing Galactic and extragalactic 
distance scales. The Preston index and spectroscopic obser- 
vations are used to determine their atmospheric metallicity 
[M] and interstellar reddening E(B — V). Since confirmed 
RR-type components are not known in binary systems, the 
mass determination is based on both stellar evolution and 
pulsation theories. Due to the uncertainty of parallax data, 
the Baade-Wesse link (BW) m ethod is mostly used to infer 
their distance d (|Smithlll995l V In the BW analysis, and to 
determine [M] and E{B — V), the quasi-static atmosphere 
approximation (QSAA) is employed to interpret photometry 

and spectroscopy. 

The QSAA was introduced by iLedoux fc WhitnevI 

l|l96Cll : 'The simplest approach is to assume that at each 
phase, the atmosphere adjusts itself practically instanta- 
neously to the radiative flux coming from the interior and 
to the effective gravity g^ 



gc = GMR^' + R 



(1) 



where R, and R are the instantaneous values of the radius 
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and acceleration, which is supposed uniform throughout the 
atmosphere', G is the Newtonian gravitational constant, M 
is the stellar mass and dot is a differentiation with respect 
to time t. 'One may then build a series of static model atmo- 
spheres,' and select one of them at each phase by spectro- 
scopic or photometric observations. Its flux, colours, effec- 
tive temperature To and surface gravity g^ are accepted as 
the atmospheric parameters of that phase providing a basis 
for the determination of other parameters such as angular 
radius, mass, distance, etc. 

The subject of the this paper is the QSAA and its 
generalization to a non-uniform atmosphere. We will inves- 
tigate the QSAA from the point of view of atmospheric 
emergent flux and hydrodynamics. By comparing the ob- 
served colour indices with those of s tatic model at mospheres 
(|Castefli. Gratton fc Kurucj|l997l . lKurucj|l997l ') we select 
the phases when they coincide. Considering hydrodynamics, 
we do not construct a consistent dynamic model of an RR 
atmosphere. However, we flnd a better description of the 
pulsating atmosphere if we characterize it by pressure and 
density stratiflcations in addition to the two parameters R 
and R. We determine the fundamental parameters Ai and 
d from the hydrodynamic considerations without using the 
BW method or theories of stellar evolution and pulsation 
at all. Our method uses photometry as observational in- 
put; spectroscopy and radial velocity observations are not 
needed. 
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In Section 2, conditions of the QSAA are formulated for 
a spherically pulsating compressible stellar atmosphere with 
a velocity gradient. Practical methods are described to de- 
termine [M], E{B-V), T^((p) and Xogg^^ip) from UBV{RI)c 
colours of Kurucz atmospheric models, ip being the pulsa- 
tion phase. The laws of mass and momentum conservation 
are used to determine the mass and distance of the pul- 
sating star. Section 3 presents the results obtained from the 
UBV{RI)c photometry of SU Dra. The discussion and con- 
clusions are given in Sections 4 and 5, respectively. 



2 THE QUASI-STATIC ATMOSPHERE 
APPROXIMATION 

First, we describe a photometric method to select the phases 
in which QSAA can be regarded as a good approximation 
from a point of view of fluxes. Next, we introduce the time- 
dependent pressure and density stratifications of the photo- 
metrically selected static models in the laws of mass and mo- 
mentum conservation to find the phases in which QSAA can 
be regarded as a more or less good approximation from a hy- 
drodynamic point of view as well. Using the phases of valid 
QSAA from photometric and hydrodynamic points of view, 
we determine mass, distance and infer the internal motions 
of the atmosphere with respect to the stellar radius R. In this 
paper, we define R as where the optical depth is zero, since 
the emergent fluxes of a theoretical model are given for r = 
in any photometric band. Of course, the zero boundary con- 
dition is a mathematical idealization (jLedoux fc WhitnevI 
1 196(3 ). In the following, we tacitly assume that the density 
is not zero at the boundary, i.e. q{R) >OifO<r<Cl. 

2.1 QSAA from the photometric point of view 

The functions n[CIi,Ch,[M],E{B - V)] and 
ge[Ch, Ch, [M], E{B - V)] of the Kurucz atmospheric 
models form a suitable grid for interpolation in the ranges 
6000 < Tc < 8000 and 1.5 < log go < 4.5 if the colour indices 
CIi,Cl2 are selected appropriately from UBV{RI)c pho- 
tometry. For fixed [M] and E{B — V) the intersection of the 
two functions {T!;'\log g„ Ch, Ch, [M], E^B - F))},=i,2 
as a function of log go gives a pair of Te(v3),logge(v) for a 
given pair of CIi((p), Cl2(¥'). Although it is more practical 
to us e U - 2B + V instead oi U - B (|Barcza fc Benkd 
l2009f ). the other common colour indices, e.g. U — V etc., 
will also be kept. 

Using UBV{RI)c photometry, we can construct four 
independent colour indices CIi=i,2,3,4, i.e. y^) ~ 6 combina- 
tions of colour indices resulting in six pairs oi Te{^p) , log ge{ip) 
for a given ip. Since the U band covers the Balmer jump, the 
indicator of gcif), the combinations must contain at least 
one colour index containing U photometry. (2) ~ (2) ~ "^^ 
such combinations can be constructed. The average of the 
30 pairs of Tc{'fi), log gc{'fi) will be accepted for a phase (p. 

Condition I. If the scatters Aro((/3), Aloggc('/') of 
Tc{p),loggc{p) are compatible with the expected scatter 
from the error of the colour indices, QSAA provides a good 
approximation in this phase from the photometric point of 
view. 

The bolometric corrections and UBV{RI)c physical 
fluxes of a model with given Tc,loggc,[M], E{B — V) are 



available in tabular form (|Kurucall997l ) and they allow us 
to determine the angular radius of the star defined by 

^(V^) = R{p)/d. (2) 

Technical details of finding the appropr iate model an d de- 
termining i9 at a phas e (p are given in iBarczal l|2003h and 
iBarcza fc Benkd l|2009t )Fl 

On the scale provided by the Kurucz atmospheric 
models, the above procedure offers a possibility to de- 
termine reddening and atmospheric metallicity from pho- 
tometry without using spectroscopy. In the shock free 
phases, the QSAA is expected to reproduce the colours 
well, therefore, ATcAlog^o must be minimal as a func- 
tion of the (/)-independe nt model parameters E{B — V), [M] 
(jBarcza fc Benkd l2009l ). Of course the method can also be 
applied for non- variable stars to determine E{B — V), [M], 
Te, logs, i9. 



2.2 QSAA from the hydrodynamic point of view 

In the original form of the QSAA (|Ledoux fc WhitnevI 1960l '). 
the atmospheric motion is described by the parameters 
R(ip),R{ip),R{p) and the outward or inward accelerations 
are driven by the variable gc{p) — gslRi^p)] in ([1} where 
(?s ~ GM/r^. This is essentially a uniform atmosphere ap- 
proximation (UAA). However, the atmosphere of an RR star 
is a dilute, compressible gaseous system with variable tem- 
peratures. Therefore, the introduction of additional param- 
eters containing Tc{ip) promises a better description in the 
frame of the QSAA. 

Conside ring the geometry of radial pulsation modes 
(|Smithlll995l ) and neglecting rotation will result in all com- 
ponents of the vector field velocity v{r,t) being zero except 
for the radial component v(r, t). The viscosity is negligible at 
the velocities occurring in an RR atmosphere. Therefore, the 
momentum conservation is expressed by the Euler equation 
of hydrodynamics: 

^" , ^"1 ^ ^ , 1 ^P I (tang) / ,N n fO\ 

where t = Pp, P is the pulsation period. 

The contribution of the neglected tangential motions is 
the average af*""'^) = (Atv)-'' J^ d O sin 9 /„ d0[u9 {dvr/de) + 
v^ sin" 



9{dvr/d(j}) -vl - vl]/r (|Landau fc Lifshitj|l98(]| '). 
In perfect spherical symmetry a'*^"^' = 0. The actual pres- 
sure and density stratifications were devided into two parts 
p(r, t) = p(r, t) -Hp('^>'"> (r, t) and ^(r, t) = Q[r, t) + f-''^^"' (r, t) , 
where p{r,t) and Q{r,t) denote the values of a static model 
atmosphere, while p^'^^^'{r,t), p''*''"' (r, f ) are the dynamical 
corrections due to pulsation. 

From the hydrodynamic point of view, the essence of 
the QSAA is that at a given t 

I dp 



g dr 



^ ffc(t) > 0, 



(4) 



i.e. the r-independent effective gravity of a static model 
atmosphere is introduced in ((S]) and p''*''"' (r, t) — 0, 
p(dyn) j-^^ ^-j _ Q j-jj ^jjgj^ follows, we drop this latter restric- 
tion and define a dynamical correction term of acceleration 

^ A program package to perform the conversion of colour indices 
and magnitudes to To, logg, i9 is available upon request. 
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(5) 



(dyn) ^ 1 ^ _ i ^ 

Q dr Q dr 

to account for the difference of a static and dynamical atmo- 
sphere. Because of the lack of dynamical model atmospheres, 
(j(dyn) jjjyg^ jjg determined empirically from the quantities 
derived from photometry. 

Equation ([3]) is transformed to definition ((l| of QSAA 
if dv/dr = and a''''^'^' = 0, a''*^") = 0, because these 
simplifications imply dv/dt = R and ga{R) = GMR~^, 
i.e. the UAA is valid. As a dynamical equation, this sim- 
plified equation of motion is the basis of mass dete rmi- 
nation in BW analyses (e.g. Cacciari, Clementini & B used 
ll989l . lLiu fc Janeslll990l .etc'). 

If the sound velocity is assumed as an upper limit for ve 
and v^, a''^"^' is some O.lms"^. It is negligible in compar- 
ison with the other components of the acceleration at any 
r '^R. Typical values are Qb < lOms"^ during the pulsation 
cycle of an RR star while Qc can exceed lOOms"^ when the 
atmosphere is in the state of maximal compression. Thus, 
(j(tang) ^jjj j^g neglected and 



g{r,t) =g4t)~gs{r) 



(6) 



is a periodic function of t with zero-points. 

If the r-dependence of the temperature T is neglected 
and a constant density upper boundary condition is as- 
sumed, an integration of (|4]) over the interval [r, R] with 
the equation of state of a perfect gas gives the approximate 
density stratification of the model atmosphere: 



gir, t) = g{R) exp{-/io(i?, t) [r - R{t)] } 



(7) 



where ho{R,t) — fj,ge{t) /TZT{R), TZ, fi are the reciprocal 
barometric scaleheight at R, the universal gas constant and 
average molecular mass, respectively. Essentially, we use the 
static model atmospheres to measure —g~^{dp/dr). The at- 
mospheric pulsation is driven by g{r, i) — a'-'^-*'"' and the ther- 
mal processes are represented by the variable ho. 

Condition II. If dv/dt + vdv/dr ~ g, 
i.e. \a'-'^^'^\r,t)\<2Ag, QSAA provides a good ap- 
proximation in this phase from the hydrodynamic point of 
view. 

This condition formulates the fact that the dynamical 
excess of acceleration in the upper photosphere is smaller 
than the Ag error of g and the same error is assumed for 
dv/dt + vdv/dr. The combination of (0) and ((3| allows us to 
estimate the constant and 0(d) terms of a''*^" . However, the 
satisfaction of Condition II can be checked afterwards when 
the terms with other powers of d were determined. The quo- 
tient q{r, t) = «''*■*'"'/(? characterizes the degree of excellence 
of the QSAA. The QSAA is exact from the hydrodynamic 
point of view ff a''*^") = 0. 

The perfect spherical sym metry means that in a Euler 
picture (|Pringle fc Kind [20071 ') we have to use ©-Q with 
(j(tang) ^^^ t^ — Q and the mass conservation law 



dg dg dv 2 
-^+VT^ + g^ + -vg^O 
at or or r 



(8) 



must be taken into account. With assumption ((Tjl, the ana- 
lytic solution to ((HJ is 



v{r,t) 



-air + ao + a_ir + a_2»' 



(9) 



where ai — Ch^^ ^{dho/dt), ao = R + Cai{R — S/Iq ^), a_i = 
2CaohQ^, a-2 = Ca-ih^^, C = I. C was introduced to 
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Figure 1. Te, i?, logge, ho{R) of SU Dra as a function of phase. 
Filled circles indicate the phases when Condition I is satisfied. 

incorporate UAA by taking a''^^'^\R, t) = and C = 0. The 
main term in u is 7? and the convergence of Q is excellent 
because rho S> 1. 

After differentiations of v, ho and i9, taking r — R, the 
following form is convenient for a numerical solution: 



M{d,ti) ~ M{d,t2) = 
where M{d,t) 



(10) 



[gdt) - (dv/dt) - v{dv/dr) 



{R,t)]R /G. It can be solved if there exist two or 



j(dyn) 

more t intervals satisfying Condition I and a^'*^"'(i?, i) 
identical for them. The assumed differences a^'^^'^'{R,ti) — 
Q,(dyn)|-^^^^'j ^ Q ^^^ ^Yie satisfaction of Condition II must 
be verified afterwards. Equation (|10p must be solved for 
different phase pairs {ipi,(f2) and the values d and M 
must be omitted from the final averaging if a^'^^'^^{R, (pi) y& 



,(dyn) 



{R,^2) 



To account for the difference of the local and effective 
temperatures, the boundary tempera ture T(R) = 2 ~°'^^Tc 
of a gray atmospheric model at r = l|Mihalaslll978h will be 
used in ho{R,t). Hydrogen and helium are mainly neutral 
at the boundary temperature of RR stars. Consequently, 
assuming a typical chemical composition means that jj, — 1.3 
is the appropriate choice in ho{R,t). 



3 RESULTS FOR SU DRA 



iBarczal l|2002f ) derived good quality light curves of the RRab 
star SU Dra by homogenizing 45 years of UBV{RI)c ob- 
servations (see his table 7). The method described in the 
previous section will now be applied. Instead of smoothing 
the light curves, binned data will be used to keep the results 
as close as possible to observed quantities. The 30 colour 
index pairs containing at least one U were obtained from 
CI,|E}° = {U-2B + V,U-V,U-Rc,U-Ic,B-V,B- 
Rc.B-Ic,V-Rc, y-Ic,Rc-Ic}. Takim [M] = -1.60, 
E(B - V) = 0.015 (|Liu fc Janeslll990[ ) the basic quantities 
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Figure 2. Upper panels: histogram of 30 log go, Te values in bins 
of size 0.1, lOK, respectively. Dotted: SU Dra at (/? = 0.5, solid 
line: SU Dra at i/? = 0.98, dashed: BD +67°708. Lower panels: 
standard errors A log g^ and ATe from 30 possible combinations 
of Cli as a function of phase. At phases below the dotted lines 
the scatter can be attributed to observational errors of the colour 
indices. 



for (|10[l are summarized in Table [T] for different phases. Fig. 
1 is a plot of Tc{(p), i?(v3), log3e(v), ho{R,ip) for one whole 
pulsation. 

ATc(tp) and Alog(;e(<^) are plotted in the lower panels 
of Fig. [2] Assuming random errors of ±0.02 for the colour 
indices will result in ATe = ±10K and Alog^e = ±0.04. 
These values are indicated by the dotted horizontal lines. 
Condition I is satisfied in the phase points lying below or 
close to the dotted lines. At phases lying above the dot- 
ted lines, the monochromatic fiux of static models and SU 
Dra differs significantly on a level which has noticeable ef- 
fect on the broad band colours UBV{RI)c- The assumed 
[M] = —1. 6, E{B — V) — 0.015 were verified by a variation 
procedure (jBarcza fc Benko 2009 ). In the shock free phases 
ifi = 0.15,0.5,0.55, minimization of Are((p), Alog(?e(v5) re- 
sulted in [M] = -1.60 ± 0.10, E{B -V) = 0.015 ± 0.01. 

To demonstrate the difference between good and poor 
QSAA, the histograms of the 30 log (jc, T^ values are plotted 
in the upper panels of Fig.[2]for </3 = 0.5, 0.98 of SU Dra and 
BD ±67 708. They show normal distributions with small 
scatter for the non- variable BD ±67° 708 and SU Dra at 
ifi = 0.5. At </5 = 0.98, the distribution is almost uniform. 
At the next phase point ip = 1, merely 14 intersections of 
{ri'^(log<?e,CIi,Cl2,[M],£;(B - V))h=i,2, i.e. 14 pairs of 
Iog(7e,Te were found instead of 30 pairs. Thus, at (p ^ 1 the 
observed colours differ significantly from those of any static 
model of iKuruca (1997). However, it is interesting to note 
that the small errors Alog^c ~ 0.03, ATe — 12K indicate 
a phase island at ip = 0.93 with observed and static model 
colours in agreement for some O.OIP « 10 minutes. This 
happens to be the phase of the hump on the light curve when 



Table 1. The basic quantities for mot at various phases. The 
units are cms~^,K,rad X 10^°,cm~^ X 10^". 



v 


loggc 




Te 


^ 




ho{R) 




0.1 


3.08 ± 


.07 


7418 ± 27 


1.627 ± 


.004 


30.0 




0.15 


2.59 ± 


.03 


7021 ± 10 


1.723 ± 


.003 


10.4 


I,II 


0.2 


2.16 ± 


.04 


6735 ± 14 


1.792 ± 


.004 


3.99 


I,- 


0.25 


2.03 ± 


.04 


6596 ± 12 


1.801 ± 


.004 


3.00 


I,II 


0.3 


2.14 ± 


.04 


6519 ± 13 


1.799 ± 


.004 


3.96 


1,11 


0.35 


2.39 ± 


.07 


6495 ± 18 


1.774 ± 


.004 


7.10 


-,n 


0.4 


2.67 ± 


.08 


6524 ± 21 


1.723 ± 


.004 


13.4 


-,- 


0.45 


2.70 ± 


.05 


6519 ± 13 


1.692 ± 


.004 


14.3 


-,- 


0.5 


2.57 ± 


.01 


6418 ± 3 


1.717 ± 


.003 


10.7 


I,- 


0.55 


2.38 ± 


.04 


6328 ± 8 


1.740 ± 


.003 


7.11 


1,11 


0.6 


2.29 ± 


.06 


6292 ±11 


1.741 ± 


.004 


5.83 


-,- 


0.65 


2.30 ± 


.06 


6374 ± 11 


1.678 ± 


.004 


5.77 


-,- 


0.7 


2.25 ± 


.06 


6448 ± 17 


1.616 ± 


.006 


5.18 


-,- 


0.92 


3.95 ± 


.06 


7055 ± 22 


1.546 ± 


.006 


235 


-,- 


0.93 


4.05 ± 


.03 


7300 ± 12 


1.545 ± 


.003 


289 


I,- 


0.94 


4.17 ± 


.05 


7502 ± 27 


1.520 ± 


.006 


371 


-,n 


0.95 


4.01 ± 


.07 


7582 ± 42 


1.520 ± 


.008 


164 


-,- 



Note. Satisfaction of Conditions I and/or II is indicated in the 
last column by I and/or II, respectively. 

the in ward and ou tward motions encounter and produce a 
shock (|Smithlll995l '). 

In the interval 0.92 < if < 1.05 the atmosphere is in 
a state of maximal compression by the shock coming from 
the sub-photospheric layers and R is nearly minimal. This 
is the rising branch and the start of the descending branch 
in the light curve. Therefore, the values of log(;e,Te,i? ob- 
tained from QSSA, if they can be found at all, must be 
considered as a first approximation only. This is reflected in 
large Alogpc, ATo except for ip ~ 0.93. 

To obtain v{r,t), etc., for (|10|) . i9(i) and ho{R,t) were 
differentiated by midpoint formulae. Fig. El^a) is a plot 
of the functions Ai{d,t) for the phases t = ipP, ip = 
0.15-0.35,0.5,0.55 with a<'*y"'(i?,f) = 0. The average and 
standard error oi M,d are given in Table [2] from the pairs 
[ifii = 0.25 and ifi2 = 0.15,0.3,0.35,0.55), (ipi = 0.35 and 
V32 = 0.55) in (I10|l as our best values denoted by [*]. At 
ip — 0.35 Condition I is moderately violated, but Condition 
II is satisfled and a'-'^y''\R,t)/gs{R,t) ^ -0.07; therefore, 
this phase was included to obtain d, M of [*] . Condition I 
is satisfied at (^ = 0.2, 0.5; however, these phases had to be 
excluded from the mass and distance determination because 
of the large a'^'^^''\R,ip = 0.2,0.5) = 5.8, -3.3 ms^^ and 
g = —0.52, 38, respectively. 

Using our best solution for M and d, the radius vari- 
ation, velocities and the components of acceleration were 
computed in physical units and are plotted in Fig.[3Ib)-(d). 
Velocities and accelerations are plotted only for the phases of 
more or less good QSAA {tp — 0.15-0.3,0.5,0.55), including 
the slightly shocked phases (p — 0.35-0.45. 

At the phases ip = 0.15,0.25,0.3,0.35,0.55 
a''*y">(i?,(/9) = -79,8, -16,65,93 cms~2 and Condition 
I is satisfied, \a'''^''''\R,t)/gs{R,t)\ < 0.13, i.e. a^'^'"'\R, ip) 
is small in comparison with the other acceleration terms 
in (IS}. |g| ~ 0.1 is expected from Alogjc = 0.04, 
q{R,'p) = 0.39, 0.02, 0.06, -0.10, -0.13 were found. 
Thus, Condition II is indeed satisfied. The outlier value of 
0.39 is produced by cancellation of dv/dt ~ —30.5 ms~^ 
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Figure 3. Panel (a); The functions M{d,ip) in units of so- 
lar mass and pc, respectively. Dotted, from left to right: ip = 
0.2,0.5, dashed: ip = 0.25, lines: from left to right ip = 
0.55,0.35,0.3,0.15. Panel (b): crosses: R{f), squares: R((p) — 
l/ho{R), in solar units. Panel (c): crosses: R{f), squares: v{ip), 
triangles: h^ [dv/dr)\j.—ii. Panel (d): crosses: R(^), squares: 
dv/dt + vdv/dr, circles: a('^^°)(i?, ip), thick line: ga{R). In panels 
(b)-(d) the angular variables were converted to absolute radius 
etc. by d = 663 pc and M = O.68X0. 



Table 2. Distance and mass of SU Dra from ljlO|l . 



d[pc] 


M[Mq] 


Used phases 


Remark 


663 ± 67 


0.68 ±.03 


0.15,0.25-0.35,0.55 


c = i, w 


658 ± 104 


0.65 ±.03 


0.15,0.25,0.30,0.55 


C = 1 


618 ± 200 


0.54 ±.10 


0.20,0.30,0.50,0.55 


C = 0, UAA 



Note. [*] indicates the values which we accept as the best ones. 

and vdv/dr ~ +28.2 ms~^, i.e. the free fall of the upper 
atmosphere is for a short time almost uniform with constant 
velocity and vanishing acceleration —g{R,(p) ~ — 2.3 ms~^. 
Therefore, the density stratification differs from (O and 
this is manifested in the large q. However, this phase could 
also be used to get the best d and M in Table [2] because 



,(dyn)| 



< 1 ms 



In the interval 0.92 < f < 0.95, the acceleration 
is dominated by g^ ~ 89, 112, 148, 102 ms~^ S> Qs and 
the atmosphere is practically at a standstill. The values 
q = 1.06,-1.68,-0.12,1.44 and a'^'^^"-\R,t) / ga{R,t) « 
5,28,-2,6 indicate a'^^y") si 41, 254, -19, 55 ms"^ and 



,(dyn)l 



'^ Qs- If <?s is neglected by setting M(d,t) — and 



the data of Table [T] are used with 



Jdyn) 



0, the upper limit 



d < 642pc is obtained from (p — 0.94. If this estimation is im- 
proved by a''*^"'(i?,v5 = 0.94) = -19ms"^ obtained for the 
distance and mass [*], the upper limit shifts to d < 724pc. 
From the functions M{724pc, <fi = 0.15,0.25-0.3,0.55) the 
upper mass limit is A^ < 0.84 ± .03 Mq. These upper lim- 
its are in accordance with the values [*]. However, they 



should be treated as indicative only, because they are sensi- 
tive to the actual value of 1?, 1?, dho/dt, d^ho/dt^. These val- 
ues could be determined by numerical differentiations from 
i9, ho in 0.91 < <^ < 0.96 when the QSAA is not optimal to 
obtain all of them. Hence, the phase island (p — 0.93-0.94 
could not be included in the mass and distance determina- 
tion because a very transient good QSAA is embedded in a 
strongly shocked interval and the mass Ml appearing in Qs 
plays a secondary role in the dynamics. 

For the sake of completeness, it must be mentioned that 
the pair {(pi — 0.15 and (p2 — 0.3) in (|10l) gives d = 1485 pc. 
Adding this value to those involved in [*], d — (800±148) pc 
and A4 — (1.15±.O5)A40 is obtained. This distance is closer 
to d = 900 pc suggested by t he Hipparcos parallax vr = 
1.11 ±1.15 mas (|PerrvmarJll997l ) . However, this 1485 pc was 
omitted because it is an outlier point above the 3cr level, it 
exceeds the upper limit 724 pc and the inaccurate knowledge 
of a('*y"'(i?,0.15) - a('*y")(i?,0.3) makes it uncertain. 

Close to the shocked phases (i.e. for ip = 0.1, 0.4, 
0.45, 0.7-0.76, 0.84, 0.86, 0.91, 0.94, 0.97), negative stel- 
lar masses were obtained if a-'^^'^' {R, ip) = 0. The condition 
M{d,ip) > gives an upper limit for a''^^"'(i?, ip) < at 
these phases and shows clearly the limitations of treating an 
RR atmosphere in the QSAA. In these phases the character- 
istic value \q\ « 1. The sharp peaks q{R, ip = 0.5) = 38.8 and 
q{R, p> = 0.65) = —13.8 are produced by dv/dt « —vdv/dr, 
i.e. the atmosphere is almost completely free of acceleration. 
q[R,ip = 0.05) = 11.4 is a result from the poor representa- 
tion of the density stratification by ([7]), i.e. the dynamical 



term is large in Qc — Qs 



,(dyn) 



= (44.7-8.7-33.1): 



producing a<'*^"'(i?, (^) ^ ll{dv/dt + vdv/dr). This is the 
start of the descending branch in the fight curve, when the 
atmosphere starts expanding rapidly. 

The averaged surface gravity for the whole pulsation 
cycle is (loggc) — 2.73, which is in excellent agreement 
with (log fle) = 2.69 2.72 from uvby an d UBV{RI)c pho- 
tometry (jSiegell 11983 . iLiu fc Janej|l990l 'l. respectively. The 
averaged effective temperature of (Te) = 6778 K is signifi- 
cantly hig her than (Tc) = 6433,6400,6490 K feiegell ll982l . 
iLiu fc Jan es 1990, Barcza 2003). The minimal, averaged and 
maximal angular radii are (1.519, 1.669, 1.801) x 10"^" rad 
corresponding to i?min = (4.46, 4.90, 5. 29)7?©, respectively. 
The follo wing equilibrium luminosity and effective tem- 
perature (|Carnev. Strom fc Jones Ill992 i) give the position 
of SU Dra in a theoretical Hertzsprung-Russell diagram 
(HRD): Leq = 47rcrd2(^2(<^)r*(v9)) = (45.9 ± 9.3)1,0, T^q = 
{Leq/47rcr[(i9)d]2}i/4 ^ ^^2y4^i/4^^^-i/2 ^ (6813±20)K. 

The magnitude averaged absolute brightness is {My) = 
+0.68 ± .23 mag. 



4 DISCUSSION 

From the point of view of radiative transfer, th e plane- 
paral lel approximation of the model atmospheres (JKuruca 
Il997i ) could be applied beyond doubt, because h^^ < 0.03J? 
holds for the whole pulsation cycle of an RR star. In favour 
of applying QSAA, semi-quantitative arguments were that 
the temperature changes are very slow even in a high am- 
plitude RR star such as SU Dra: 5K/2500s < \dTe/dt\ < 
1500K/5000s, while the charact eristic times-cale of the ra- 
diative processes is below 200s (|Oke. Giver fc Searlelll962l . 
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iBuonaura et al.lll985h . The effect of the variable effective 
gravity can be characterized by our Condition II, which pro- 
vides information on whether QSAA may be assumed. The 
hmits are 4.5cms-V2500s < \dge/dt\ < 14200cms-2/2500s. 
The upper hmit comes from the rising branch of the hght 
curve, when the main shock hits the atmosphere. 

In general, the photometric input of the present method 
is identical with that of the BW method. In order to obtain 
the fundamental parameters, radial velocity data and their 
problematic conversion to pulsation velocities are not nec- 
essary. On the other hand, -d and ho must be differentiated 
numerically; differential quotients are sensitive to the non- 
validity of QSAA. 

Our photometric and hydrodynamic considerations re- 
vealed empirical quantitative conditions to find the phase 
int ervals of the pulsation when the static model atmospheres 
of iKurucj 1 1997) are satisfactory to derive the variable 
and non-variable physical parameters of the pulsating at- 
mosphere. Outside these intervals, dynamical model atmo- 
spheres are necessary to refine the parameters from QSAA, 
which is beyond the scope of this paper. The fundamen- 
tal parameters d, M, [M], E{B — V) were determined using 
photometric quantities only in phases when the QSAA was a 
good approximation (i.e. both Conditions I and II were sat- 
isfied). The values obtained from averaging over the entire 
pulsation cycle can be considered as a first approximation 
only, because QSAA was assumed in all phases regardless of 
it being a good or poor approximation. The large error of 
Leq and (My) originates from Ad/d ~ .1 of our best value 
[*] in Table il 

To give a good impression of the accuracy of invert- 
ing the UBV{RI)c photometry to physical parameters, the 
comparison star BD -1-67 708 was used, bec ause its colour s 
are similar to those of SU Dra (table 3 in lBarczall2002h . 
The results are [M] = -0.77 ± .03, E{B - V) = .000, ■& = 
(1.913±.002)xl0-^° rad, log3 = 3.59±.01,Tc = 7505±5 K. 
The errors are roughly in the same order of magnitude as 
those of SU Dra when Conditions I and II are satisfied. 



4.1 Kinematics of the atmosphere 

Figs.lHc), (d) demonstrate that significant corrections must 
be added to R, R if the true pulsation velocity and ac- 
celeration are required at < r < 1. Moreover, the tri- 
angles show another correction: — 4.6kms~^ < v{R,ip) — 
hg^dv/dr < 8.3kms-^ ff 0.1 < lys < 0.6, i.e. v{r,<p) = 
v{R, ifi) — v{R, ip)ho{R, ifi){R — r) -I- • ■ ■ is the velocity profile. 
In other words, even in the shock free phases, considerable 
phase-dependent velocity and acceleration gradients exist in 
the layers R — hg^ < r < R oi the line formation. The prob- 
lem of -u 7^ is present in all phases, i.e. even in the shock 
free intervals, which are used in modern BW analyses (e.g. 
iLiu fc Janes 1990). 

The non-uniform motion of the outermost layers intro- 
duces uncertainty when pulsation velocities are determined, 
because the centre-of-mass velocity v-, must be subtracted 
from the observed radial velocities. A rece nt exposition of 
the p roblem for Cepheid stars is given in I Nardetto et al.l 
l|2009l) . By definition i9,i9,i? refer to < r < 1, while the 
variable component of the radial velocity is an average of 
velocities (e.g. (|9])) over the layers R — h^^ ^r < R i.e. < 
r < 0.3. The effect on v-y and d has not yet been studied at 



all. Nevertheless, the importance is obvious, sinc e an error 
of 1 k ms~^ in v-y results in an error Ad/d ~ 0.1 IjGautschvl 
[19871) . There is a considerable uncertainty of v-., in the litera- 
ture, suggesting an error of d as large as a factor of 2. It suf- 
fices to ment ion that for SU Dra Vy = — 161, — 166.9kms~^ 
are gi ven by lOke. Giver fc Searld (jl962h and I Liu fc Janed 
(|l99Cl ) from high dispersion spectra and spectral masking 
method CORAVEL-technique, respectively. 

A wavy fine structure in the variation of R is clearly 
seen in Figs. [T] [3jb). The outward motion starts at <,o ~ 
0.45, 0.74, 0.94-0.98 while the outward motion is reversed 
ip ~ 0.25, 0.55, 0.78, i.e. the atmospheric velocity is v{r, t) ~ 
at these phases. The wavy fine structure is without doubt 
real, because the beginnings of the outward motion are con- 
nected with significant maxima in To (<^), log cie( f)- The well - 
known bump and hump at <^ « 0.74, 0.94-0.98 (ISmithlll995h 
in the light curve are caused by the precursor and main 
shocks, respectively. At ip fa 0.45, a small change is observ - 
able in the slope of the light curve (JBarcza fc Benkoll2009l ). 
Here, the existence of a shock is a new finding. It is a pre- 
precursor shock; we propose the designation jump for it. By 
integrating radial velocities, authors ten d to smooth out the 
mentioned fine structure of motions [e.g. lLiu fc Janea (119901 ) 
in the case of SW And]. This practice seems to be unjusti- 
fied. 

The velocity of the atmosphere is v{r,t) ~ in the 
interval ip « 0.82-0.90, the atmosphere is roughly at a 
standstill, the brightness starts rising, a small depression 
of 'i}{(p),\ogge{ip) is visible a,t ip Ki 0.84-0.86, and Tg(<(S) is 
monotonic. It is not clear whether the small undulation of 
t?(<^) around (1.571-1.579) x 10"^" is real or not, because the 
maximal Alog gc ~ 0.116 was found just at <p = 0.88. 



4.2 To scale, mass, distance, absolute brightness 

Surprisingly, in comparison with previous studies, consider- 
able differences were found only in (Te), Toq. This is due to 
the fact that the interpolated Tf.{<fi) depends on log(je(v5) 
and that the information from a five colour photometry was 
used in a more complex manner: an averaged value from 30 
colour index pairs was utilized. The photometry covers the 
whole spectrum between 350 and 1000 nm. The use of only 
one colo ur index, e.g . V — K solely 'because of its apparent 
merits' (JLiu fc Janes 1990), can result in a systematic error 
of Tcif). It can explain the ~ 350 K difference, because 
in that previous work QSAA was assumed in all phases 
despite violating both Conditions I and II. An inspection 
of the functions jT^''' (lo g fle, CIi CI2, \M], E(B - V))}^^l,2 
derived from the tables (|Kurucall997l ) shows that if merely 
one colour index is to be used, the optimal choice for deter- 
mining Tc would be Re — Ic- This is because Re — Ic is 
almost independent of log g^ for the actual values of [AI] and 
E{B — V) of SU Dra. However, in phases violating Condition 
I, the sole use of Re — Ic would also introduce a systematic 
error, similar to the use oi V — K. 

There is a remarkable decrease of Ad = 200 — >■ 67pc, 
AM = 0.10 -^ O.OSMq if our compressible QSAA is sub- 
stituted for UAA, i.e. more physical input is used in the 
frame of a one-dimensional model in space. Since (|10|) was 
derived from the dynamic equation ([HJ, the mass determi- 
nation is more accurate from this while the distance is more 
uncertain. This is reflected in the shape of the curves in 
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Fig. |3ja). The large formal error of d originates from the 
non-separability of errors in i?, ■&, •&, ho, dho/dt, d^ho/dt^ and 



,(dyn) 



(r,t). 



A first attempt to derive d and A4 of an RR st ar from 
the ov er-simplified version ([l]) of ((3| was described bv lBarczal 
l|2003t) . The present results for the distance and mass of SU 
Dra differ only slightly from the previous d = (647 ± 16)pc, 
A4 = (O.66±.O3)A^0. The most probable reason of the very 
good coincidence is the existence of the phase island with 
the good QSAA at (^ ~ 0.93, just when the atmosphere is 
at a standstill; consequently, ([l]) is a good approximation 
because of v « 0. 

The values d = 640 pc, M = 0A7Mq given by 
ILju fc Janed (|l990l ) are very close to the present ones. How- 
ever, some caution is appropriate because of the underesti- 
mated uncertainties in their derivation. Their Ml originates 
from using U]), i.e. C = 0, the UAA corresponding to line 3 
in Table (2] Furthermore, because of the uncertain value of 
v-y, Ad/d<0.6 is well possible. 

In the BW method, the propagation of the error Av^ = 
1 kms~^ can be estimated from 1^ as follows. Typical ra- 
dius changes of an RR star are -Rmax — ^min ~ 5 x 10^ km 
within P/2 ~ half day. The error of the radius change is 
AR f» Av^P/2 f» 5 X 10" km, i.e. A7?/(i?max - -Rmin) « 0.1. 
By d = (-Rmax - -Rmin ± Ai?)/(^niax - i?min), the final er- 
ror will be Ad/d ~ 0.1. This considerably ex ceeds the er- 
ror o riginating from the projection factor (e.g. iLiu fc Janed 
[l990|) converting the observed radial velocity to pulsation 
velocity. Furthermore, by © and (O, the error Ad/d « 0.1 
propagates to an error AM/Ai « 0.2. 

The lower limit of the magnitude averaged visual ab- 
solute brightness is AIv > -1-0.48 mag if d < 724 pc, while 
Mv = -1-0.26, -f 0.01 mag belong to the improbable values 
d = 800 and 900 pc, respectively. 



5 CONCLUSIONS 

Observed colours and magnitudes of a spherically pulsat- 
ing star have been compared with those of static Kurucz 
model atmospheres to determine fundamental parameters 
of the star in the frame of the QSAA. Photometric and hy- 
drodynamic conditions have been formulated for the validity 
QSAA in spherically pulsating stars. 

(1) The QSAA has been generalized for a non-uniform, 
compressible atmosphere with radial velocity gradient. This 
is a step forward because the hitherto available UAA de- 
scribed the motion of the atmosphere by effective gravity 
and differentiating the sole parameter radius R of optical 
depth zero with respect to time. 

(2) Our combined photometric and hydrodynamic method 
uses the variation of effective gravity, angular radius, veloc- 
ity, acceleration in the Euler equation. The complete input 
of the Euler equation has been derived from photometry and 
theoretical model atmospheres. Spectroscopic and radial ve- 
locity observations were used neither explicitly nor implic- 
itly. This is a definite advantage in comparison with the BW 
method, because less observational efforts are needed to get 
the fundamental parameters and the problematic conver- 
sion of observed radial velocities to pulsation velocities is 
not necessary. 



(3) Concerning the interpretation of multicolour photome- 
try, the inputs are identical with those of the BW method. 
The present refinements are the quantitative conditions 
whether photometry can or cannot be interpreted by static 
model atmospheres. 

(4) First, the phases have been selected by the photometric 
conditions when the QSAA is valid. Secondly, in a number of 
these phases, the laws of mass and momentum conservation 
have been applied in the Euler formalism of hydrodynam- 
ics to determine the mass and distance of the star from the 
motion of the atmospheric layers in the neighbourhood of 
zero optical depth. Afterwards, it was checked whether the 
hydrodynamic condition of the QSAA was satisfied in the 
used phases; i.e. phases that satisfied the photometric con- 
dition but violated the hydrodynamic condition had to be 
excluded. Atmospheric dynamical mass seems to be an ap- 
propriate term to indicate that the mass has been derived by 
a method which is completely diffe rent from a mass derived 
by pulsation or evolution theories (jSmithll 19951 ). 

(5) As a by-product, a variation procedure has been given 
for estimating atmospheric metallicity and interstellar red- 
dening toward a star from UBV{RI)c photometry. This 
method was successfully applied for the non-variable com- 
parison star BD -1-67 708, giving [M] = -0.77 ± .03, E{B - 
V) = .00, ■& = (1.913 ± .002) X 10"^" rad, log g = 3.59 ± .01, 
Te = 7505 ± 5 K. 

(6) Using the UBV{RI)c photometry of the high ampli- 
tude RRab star SU Dra, the following fundamental param- 
eters have been found: 

[M] = -1.60 ±.10, 

E{B-V) = 0.015 ±.010, 



d = (663 ± 67)pc, 



Rn 



4.46i?0 , R, 



5.29R, 



■@, 



M = (O.68±.O3)M0, 

Leq = (45.9±9.3)L0, 

Toq = (6813 ± 20)K. 

Lcq and Toq are approximate values, since they originate 
from averaging over the whole pulsation cycle containing 
phases in which the QSAA is merely a first approximation. 
(7) The internal motions of the atmosphere with respect to 
zero optical depth have been found to be significant in com- 
parison with velocities and accelerations derived from the 
UAA. From the internal motions of the atmosphere, some 
constraints have been sketched for converting observed ra- 
dial velocities to pulsation and centre of mass velocities. Es- 
timations have been given for the error propagation in a 
Baade-Wesselink analysis. 
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